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Abstract 

It is argued that a spontaneous vortex phase may exist in the recently discov- 
ered compound ErNi2L>2C at temperature below 2.3K. The consequences of 
this proposal are discussed. In particular the magnetic response of the system 

are studied both above and below 2.3K and further experiments proposed. 
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Many years ago it was proposed |1]-§J that exotic phases with co-existence of super- 
conductivity and magnetism may occur in systems with competing superconducting and 
ferromagnetic components. The analysis was based on the Free energy functional 0] 

^/,M>l 2 + >l 4 + 1^-^ + 1 

+ia|M| 2 + \l)\M\ l + \l 2 \VM\ 2 - B.M], (1) 

where B = V x A, M is magnetization and ip is the superconducting order parameter. 
It was shown that a stable spiral phase where supeconductivity co-exists with spiraling 
magnetization or a spontaneous vortex phase where magnetization is more or less uniform 
in the system but vortices are generated without an external magnetic field may occur. 
Subsequently, the spiral phase was discovered in ErRh^B^ and HoMoqS$ compounds ]5]|| 
in a narrow temperature region between a superconducting phase and a ferromagnetic phase. 

More recently, it was discovered that competition between superconductivity and fer- 
romagnetism may occur in a new material ErNi 2 B 2 C . We shall show in this paper that 
ErN i 2 B 2 C is a good candidate for the spontaneous vortex phase, or that the spontaneous 
vortex phase will become stable under a relatively weak external magnetic field. (We define 
the spontaneous vortex phase in the presence of magnetic field as a state where the density 
of vortices present in the superconductor is larger than that given by the external field) Con- 
sequences of our proposal will be studied. To begin with, we first review some basic features 
of the Ginsburg-Landau free energy functional (|l|) where stability criteria associated with 
various plausible phases is examined. The analysis of Ref.(l) and (2) are then extended to 
include the effects of an external magnetic field. 

The competition between magnetism and superconductivity appears in Eq. ([[]) as a 
Meissner effect of the superconducting component towards the internal magnetic field pro- 
duced by the magnetic component B = AnM. The existence of the spiral and spontaneous 
vortex phases in F are direct consequences of Meissner effect where a uniform magnetization 
M cannot co-exist with a uniform superconducting order parameter ip. For systems with 
superconducting transition temperature T c higher than the magnetic transition tempera- 
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ture T m , a spiral phase may be stable at temperature T s less than but around T m . The 
wavevector of the spiral is of order Q ~ (A £j\/) ^ 2 where A is the penetrating depth of the 
superconducting component and £m ~ l 2 /& is the coherence length of the magnetic part 
0. The Meissner effect is avoided by having a magnetization whose average is zero on a 
length scale much smaller than Ao- At lower temperature, a ferromagnetic state with super- 
conductivity completely destroyed is usually lower in energy because of the higher energy 
gain associate with magnetization (~ ksT m ) compare with the energy gained by supercon- 
ductivity (~ (kBT c ) 2 /Ep). Alternatively, a spontaneous vortex state where magnetization 
is uniform and superconducting component exists in the form of a vortex state may be more 
stable than ferromagnetic state because of the gain in energy from the superconducting 
component. However, this state exists only if the internal magnetic field generated by the 
magnetization B satisfies the inequality 

H cl <B ~ 4vrM < H c2 , (2) 

where H c i and H c2 are the lower and upper critical magnetic fields associated with ip. The 
inequality expresses the fact that a nearly uniform magnetic field can be sustained in a 
superconductor only when the density of vortices is such that the average distance between 
them I satisfies the inequality £o < I < Ao, where £o = \Jh 2 /2m\a\ is the coherence length of 
the superconducting component. The spontaneous vortex state is favored only in systems 
where the saturated magnetization is not too strong or too weak compared with H c /47r. 
Magnetic anisotropy also plays strong role in deciding the relative stability of various states. 
In particular, easy-axis anisotropy always disfavors spiral (or linear polarized) states over 
ferromagnetic or spontaneous vortex states. 

Experimentally it is found that ErNi 2 B 2 C is superconducting below 10.5K J7| and or- 
ders antiferromagnetically with a fundamental incommensurate wave vector of (0.553a*, 0, 0) 
below 6.0K ||. The magnetic moments reside mostly on the Er 3+ ions which has a mea- 
sured magnitude ~ 8//_b 0. M Versus H measurements indicate that the compound is 
magnetically strongly anisotropic with the Er magnetic moments essentially along only the 
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in plane easy axis in (100) and (010) directions The same measurements with applied 
field along either of the in-plane axes indicate a series of meta-magnetic (field-induced) tran- 
sitions as function of magnetic field at temperature around 2.3A" [0,0. In particular, once 
the external magnetic field is significantly larger than H c \ ~ 500G, it is found that the sys- 
tem has a ferromagnetic component. The extrapolation of M(H) data back to zero applied 
field gives a ferromagnetic ordered moment of roughly 0.33/z#/ Er ||. Zero- field specific 
heat measurement shows also a break in the slope of C Vs. T curve at T ~ 2. 3 A |7|||. 
The existence of ferromagnetic component in the system is further supported by studies 
on similar compound TbN 12B2C which does not manifest superconductivity at T > 2A, 
but has a phase transition from antiferromagnetic order to an ordered state with a weakly 



ferromagnetic component of 0.5/ie/T6 for T < 8A |H|. The magnetization versus applied 
field behaviour for the two compounds are found to be very similar except at low field re- 
gion where M is always positive in the Tb compound (no Meissner effect). The similarity 
between the two compounds suggests that the magnetic behaviour of the two compounds 
are essentially identical, in particular a similar transition to an ordered state with weakly 
ferromagnetic component is also perferred in the Er compound at around 2.3A - except 
that such a transition is forbidden by the superconducting component. The origin of the 
incommensurate antiferromagnetic and the transition to the weakly ferromagnetic transition 
is not quite clear and probably involves both the exchange interactions and the dipolar in- 
teractions which are of similar magnitude. Because of the strong easy-axis anisotropy which 
forbids smooth deviations from antiferromagnetic state, the incommensurate antiferromag- 
netic state probably consists of spin domain walls separating antiferromagnetic domains, 
and the weakly ferromagnetic state is formed by ordering of the domain walls. The distance 
between domain walls estimated from the incommensurate wave vector is around 19a, giv- 
ing effective magnetic moment of around 0.42/zg/ Er, which is close to the experimentally 
zero-field extrapolated value of 0.33/^b/ Er. 

At distance scale >> lattice spacing, the antiferromagnetic component plays a negligible 
role and the competition between superconductivity and weak-ferromagnetism can be de- 



scribed by a Ginsburg-Landau functional similar to Eq. (|l|), except that the M 4 term must 
be modifed to account for the strong easy-axis anisotropy in this material. The internal 
magnetic field created by magnetic moment of 0.33fiB/Er is approximately 500G ~ H&, 
which is found to be marginal for supporting a spontaneous vortex state. However, a rel- 
atively weak magnetic field ~ H c i should be enough to drive the system from the spiral 
state into the spontaneous vortex state. In the following we shall investigate this scenerio 
using the GL functional ([!]). We shall assume that the magnetization M lies only on the 
x-y plane and shall consider external fields only in in-plane directions. The anisotropy in 
in-plane directions is not included in our analysis. We shall consider T c > T m and shall 
concentrate on the behavior of the system around T ~ T s which is the regime of experi- 
mental interest. The possibility of the system making a second order phase transition to 
spiral state at T = T s but driven into spontaneous vortex by an external magnetic field will 
be studied. We shall also discuss the alternative possibility of the system making a direct 
first order transition into the spontaneous vortex state from superconducting state. First 
we consider the temperature region T > T m and study changes in behavior of the system 
as T — > T m . In this temperature range M is small and we can neglect the M 4 term in the 
GL functional. The qualitative behaviour of the system at this temperature range can be 
most easily understood by considering the London limit where ip = constant and neglecting 
the |VM| 2 term in F. It is then easy to minimize F with respect to M and A to obtain 
M = B/a, and A = Aq(1 — 47r/a)V x B. Putting M and A back into F, we obtain 

F ~ J d 3 r 

where Ag = mc 2 /87re 2 \if}\ 2 is the London penetration depth for the 'pure' superconducting 
component, a is a decreasing function of temperature and the magnetic transition (in the 
absence of superconducting component) occurs at a(T) = 4tt. Notice that the presence of 
magnetic component reduces the overall cost in magnetic energy of the 'pure' superconductor 
by a factor (1 — An /a). It also reduces the London penetration depth from A to A = 
(\Jl — 47r/a;)Ao. The reduction in penetration depth implies that the effective superfluid 
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density observed in experiment will increase rapidly as T — > T m . As a result the critical field 



H c goes down by the same factor Jl — An /a. An interesting consequence of the free energy 
(0) is that the lower critical field H c i is not much affected by the presence of the magnetic 
component though the penetrating depth A is strongly reduced as T — > T m . To see that we 
consider the superconducting component in the extreme type II limit A >> £o- in this limit 
the energy of creating a vortex line per unit length e can be computed using the free energy in 
the London limit (H). For usual superconductors this quantity is given in the London limit 
by eo ~ (<E>o/47rAo) 2 Zra(Ao/£o) where $o is the magnetic flux quantum, fll] In the present 
case Ao — > A which enhances e. However, e is also reduced by the same factor (1 — An/a) 
because of the overall reduction in free energy (|3|). As a result, e ~ (<3>o/47rAo) 2 ^n(A/£o) in 
present case and is only reduced through reduction in A in the logarithm. Consequently, H c i 
is also not much reduced from its clean superconductor value in the presence of magnetic 
component. It is interesting to ask what happens in the system when external magnetic field 
of order ~ few H c i is applied. For usual superconductors the density of vortices is of order 
(27rAo) _1 when magnetic field is of order ~ few H&. The density of vortices is consistent 
with magnitude of external field which supplies magnetic flux of order ~ few $o i n area 
27tAo- In the present case, for external field of order ~ few H cl , the external flux supplied in 
area 2ir\ 2 is of order H c \ x 2tt\ 2 ~ $ x (A/Ao) 2 , which is much less than one flux quantum 
if A << Ao, and seems to imply that the density of vortices is much less than (27rA 2 ) -1 in 
this case. This conclusion is in fact incorrect because the total magnetic field 'sees' by the 
superconductor B = H + AnM is much larger than H as Ao >> A. It is easy to see that 
B = H /(I — An j a) in our approximation, and the total magnetic flux the superconductor 
sees in area 27rA 2 is of order ~ few $0) implying that the density of vortices is of order 
(27rA 2 ) -1 , as in the case of usual superconductors. 

Similar analysis as above can be made in momentum space when the |VM| 2 term is 
included in the GL functional. We find that qualitative behaviour of the system is not 
modified. However, the divergence in A -1 as T — > T m is removed once the |VM| 2 term 
is included. In particular, the London penetration depth is saturated at value of order 
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A ~ (Ao£ m ) 1//2 as T — > T s . At temperatures T very close to T s , the magnetic response of the 
system is dominated by the spiral instability. We find that spiraling magnetization developes 
around single vortex solution as T — > T s , with magnitude of spiral decreasing exponentially 
as distance away from vortex core. The decaying length of spiral magnetization goes to 
infinity as T — > T s , signaling the onset of spiral instability. We find also that the energy of 
vortex line e remains finite and is given in the London limit by 

at precisely the spiral instability point where A = (Xo^m) 1 ^ 2 - Notice that A -1 is of the same 
order of magnitude as the spiral wave vector around the spiral instability. 

The behaviour of vortices at T < T s can also be studied in the London limit. In the limit 
T — > T s , the magnitude of the spiraling magnetization is small and its effect on vortices 
can be estimated perturbatively. We find that the single vortex solution is very similar 
to the solution above the critical temperature, except that the decaying length A s for the 
'extra' spiraling magnetization around vortex decreases again as T decreases below T s , until 
A s ~ (Ao£a/) 1//2 , where the perturbative solution becomes unreliable. In particular, the 
energy e for single vortex line remains of order $g/ (27rAo) through out the whole temperature 
range, with no discontinuity across the spiral transition point. 

Next we consider the situation of finite density of vortices and estimate the magnetization 
as a function of the external magnetic field. Consider the Gibb's energy functional 

G = F ~ I * T? §< (4) 

where the total magnetic field B is obtained by minimizing G with respect to B and M = 
(B — H)/4n. First consider the regime T > T m and the London limit using Eq. (|3|). Let 
the applied field strength of order ~ few H c \. The total magnetic field B can be obtained 
easily by comparing the present expression for Gibb's energy with Gibb's energy for usual 
superconductors |fTTf| . We obtain 

H^^B^mm, (5) 

a m(A/Co) 



where H Q c2 ~ $ /(2^ 2 ) and M ~ B/a - H cl [ln(H^/B)/ln(X/^ )]/ATr. Notice that the 
response of magnetic component M to magnetic field B is almost identical to response of the 
'pure' magnetic system to B except the correction term ~ —H& coming from Meissner effect. 
In particular, for small enough value of a, the magnetization may become positive in this 
range of magnetic field. The value of H c2 where superconductivity is completely destroyed 
can also be obtained easily by equating B = if c2 /(l — An/a) ~ H® 2 , obtaining H c2 ~ 
(1 — 4ir/a)H® 2 , indicating that H c2 is reduced in the presence of the magnetic component. 

The magnetization curve in temperature range T > T m thus has the following qualita- 
tive feature: (1)M = —H/Ait for H < H ci , where the value of H c \ does not depend too 
strongly on temperature, in particular there is no singular behaviour around the spiral tran- 
sition. (2)Magnetization starts to increase at H ~ H cl . For H ~ few H cl , the magnetic 
component already responds to the external magnetic field more or less as if there is no 
superconducting component in the system. (3) the magnetization continues to increase un- 
til at H = H c2 where superconductivity is destroyed. Notice that M may become positive 
already at magnetic field strength H ~ few H cl . 

Next consider the magnetization curve in the spiral phase. We find that the magneti- 
zation curves above and below the spiral transition are qualitatively similar. In particular, 
the spiral state disappears and is replaced by the spontaneous vortex phase in external 
magnetic field of the order of several times H cl . The argument is based on the observation 
that at regime of temperature T < T s , the vortex solution is not much affected by spiral 
component. In particular, the value of H c i stays more or less the same above and below the 
spiral transition temperature. As external magnetic field is of order ~ few if cl , the distance 
between vortices will be of order ~ A ~ (Ao^m) 1 ^ 2 - However, this is of the same order as 
the period of the spiral state. When the vortex distance is comparable with period of spiral, 
the spiral state losses its meaning. Thus we expect that at this magnetic field range, the 
spiral state will smoothly crossover to the spontaneous vortex state where the magnetic 
component of the system responses to external magnetic field more or less independent of 
the superconducting component as in the high temperature phase. In particular, the mag- 
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netic response will be similar to that of a 'pure' ferromagnet at temperature T < T m , in 
agreement with what is observed in ErNi 2 B 2 C compound. The properties of the sponta- 
neous vortex phase can be studied by writing M = Mo + M', where M5 2 ~ (An — a)/ (3 is the 
spontaneous magnetization of the 'pure' ferromagnetic component at temperature T < T m 
and the GL functional can be expanded to second order in M' . Neglecting the |VM| 2 term 
as before, we find M' = (B — An Mq) j \12n — 2a), and the effective GL functional in the 
London limit in terms of B and H fields has the same form as (||), except that the total 
magnetic field B is coupled to an effective external magnetic field H e ff = AnMo + r]H , where 
7] = ((6n — a) /(An — a)). The effective London penetration depth is A = ^ _1 Ao and the 
criteria for a stable spontaneous vortex state is given by 

H c2 > H eff = An M + V H » H cl . (6) 

Notice that A increases again as temperature lowers. At very low temperature, A — * Ao- 
In fig.l we show the ratio of total magnetic field to external field B/H as a function of 
external field H for H > if cl at several different values of temperatures computed using 
Eq. (^]), with corresponding equation for T < T m . We have choosen H c2 /H c i = 25.0, 
6v//£o — 0-5, Oi{T) = An + 60tt(T — T m ) and with saturated magnetization AnM = 2.0H c i at 
zero temperature in generating the figure. It is clear that the zero field extrapolation of the 
curve at T = 0.5T m indicates existence of ferromagnetic component in the system. Notice 
that B/H measures the total density of vortices in the system and is > 1 for T < T m . The 
density of vortices for the 'pure' superconductor is close to the curve with T = 3.0T m . The 
difference arriving from the ferromagnetic component is huge at T < T m , as can be seen 
from the figure. 

Lastly we want to make a few comments on the properties of the spontaneous vortex 
phase, in particular in the limit when the saturated magnetic moment is large enough and 
magnetic anisotropy is strong enough so that a direct first order transition from supercon- 
ducting phase into spontaneous vortex phase occurs in the absence of external magnetic 
field. In this case, the effective magnetic field the superconductor sees is always larger than 
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H c i and there will be no Meissner effect associated with additional external magnetic field 
applied on the system, i.e. the effective H c i of the system is zero and superconductivity 
'appears' only when vortices are pinned to impurity sites in the system. Notice that Meiss- 
ner effect exists in the spiral phase where H c \ > 0. Thus measurement of the Meissner 
effect (for example, by SQUID) will distinguish the spiral and spontaneous vortex phase 
unambigiously. Experimentally, it seems that Meissner effect are observed in the ErN 12B2C 
compound in the M vs H experiment. However, the experiment is performed in zero-field 
cooled environment |7]U indicating that the result may not reflect the true equilibrium ther- 
modynamic state of the system. Thus the possibility of a zero-field spontaneous vortex 
phase existing in the compound can not be ruled out. Direct observation of the spontaneous 
vortex phase by imaging techniques is suggested. 

We wish to acknowledge very useful discussions with Peter Gammel and U. Yaron. 
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FIGURES 

FIG. 1. B/H vs. H/H c i at H > H c i for several values of temperatures computed approxi- 
mately in GL theory. The contribution from ferromagnetic component is clear from the figure. 
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